Quantum phases of Bose-Bose mixtures on a triangular lattice 
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We investigate the zero temperature quantum phases of a Bose-Bose mixture on a triangular 
lattice using Bosonic Dynamical Mean Field Theory (BDMFT). We consider the case of total filling 
one where geometric frustration arises for asymmetric hopping. We map out a rich ground state 
phase diagram including xy-ferromagnetic, spin-density wave, superfluid, and supersolid phases. In 
particular, we identify a stripe spin-density wave phase for highly asymmetric hopping. On top of 
the spin-density wave, we find that the system generically shows weak charge (particle) density wave 
order. 

PACS numbers: 67.85.Hj, 67.60.Bc, 75.10.Jm, 67.85. Fg 



I. INTRODUCTION 

Geometric frustration arises when magnetic interac- 
tions between different spins on a lattice are incompati- 
ble with the underlying crystal geometry. Since the first 
investigation of the Ising antiferromagnet on the triangu- 
lar lattice [l|, geometric frustration has been a constant 
source of surprises that inspired the development of new 
concepts 0, yjj- This provides strong motivation to study 
the physics of frustration from the different perspective 
provided by systems of ultra-cold atoms. 

Recent experiments have made substantial strides in 
this direction with the realization of non-bipartite trian- 
gular 0, Q and kagome @ optical lattices. In this paper 
we investigate the consequence of bringing a two com- 
ponent bosonic mixture into the Mott insulating regime 
on the triangular lattice. Deep in the Mott state, the 
magnetic exchange interactions can include as the main 
source of frustration a strong Ising antiferromagnetic ex- 
change. 

Our goal is to study how the magnetic phases evolve 
when we increase the frustrating Ising interaction or ap- 
proach closer to the transition to the superfluid phase, 
such that the effective spin model with second-order ex- 
change interactions no longer applies. 

To answer this question we apply Bosonic Dynami- 
cal Mean Field Theory (BDMFT) 041, which is non- 
perturbative in the hopping amplitudes. We find that 
even in the deep Mott regime, the standard spin exchange 
model is invalid for extremely asymmetric hopping. In- 
stead, we numerically identify a stripe spin-density wave 
(SDW) phase (see Fig. QJi.2). This can be understood 
within a higher-order effective spin-model ((HJ derived 
from fourth order perturbation theory. The effective de- 
scription shows that the stripe SDW is favored by the 
higher-order density fluctuations of the "lighter" atoms 
which remove the Ising-type frustration. Moreover, we 
also find that on top of the spin-density wave, due to 
asymmetry of the hopping amplitudes, the system de- 
velops a weak charge-density wave in the total particle 
density (see Fig. [3]) . 

The paper is organized as follows. In Sec. II, we intro- 



duce the system and model studied here, as well as the 
theoretical approach used in our investigation. In Sec. 
Ill, the main part of this paper, we present a detailed 
discussion of the ground state properties of the system. 
We conclude in Sec. IV. 



II. MODEL AND METHOD 

We consider two species (hyperfine states or isotopes) 
of ultracold bosons loaded into a triangular optical lat- 
tice. For sufficiently low filling this system can be de- 
scribed by a two-component Bose-Hubbard model in the 
lowest band approximation: 

H = - y^ftqajaj + tbb\bj + h.c.) + U y^n a jnu 

(hi) » 
+ - ^ V a nia(riia - 1) - X] ^aTlai- (1) 
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Here denotes nearest-neighbor sites, aj(aj), bi(b\) 
are bosonic annihilation (creation) operators of the two 
species on site i in the Wannier representation, and n a i = 
a\a,i, ribi = b\bi. The first term in Eq. (fTJ) describes the 
kinetic energy of each species with hopping amplitudes 
t a and tb] the second and the third term represent the 
on-site inter-species interaction U and the intra-species 
interactions V a and VJ, for species a and b, respectively; 
finally, fj, a and fib denote the chemical potentials. 

Previous studies on two-component ultracold bosons 
in a square or a cubic optical lattice have revealed a rich 
phase diagram [10]. Within the Mott insulator, at low 
temperature quantum magnetism arises, in particular z- 
antiferromagnetic and xy-ferromagnetic order. For total 
filling one, the emergence of magnetic order can be eas- 
ily understood if we note that in the deep Mott regime 
(strong coupling limit) t a _b <C U, V a ,bi the physics of the 
two-component Bose-Hubbard Hamiltonian (fTJ) is given 
by an effective spin-1/2 XX Z model 
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where S, = (a\, b\){<r/2) { J with a x , a y ,a z being the 
Pauli matrices and 



a, 
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In the following discussion, we assume h — 0, i.e. van- 
ishing spin imbalance (in the case of asymmetric hop- 
ping t a ^ the chemical potentials \x a and [ib are 
tuned to achieve this). On bipartite lattices (e.g. square, 
cubic) the system supports xy-ferromagnetic order for 
J± > Jz > 0, which is characterized by the local correla- 
tor (a'b), and z-antiferromagnetic order for J z > J± > 0, 
characterized by the order parameter A a f = \(S^) — (S^)\ 
with a (a = —a) being the sublattice index. 

We expect new interesting physics to emerge in a sys- 
tem with triangular instead of bipartite optical lattice. 
For antiferromagnetic exchange coupling J z , it is impos- 
sible to minimize the energy of the spin configuration 
on each lattice bond, i.e. geometric frustration arises 
and the system may develop exotic phases at low tem- 
perature. We are particularly interested in the extremely 
asymmetric hopping regime for the full range of couplings 
from Mott insulator to superfiuid. Since higher order 
density fluctuations of the "lighter" species are expected 
to become even larger than the low-order density fluctua- 
tion of the "heavier" ones, their interplay with Ising-type 
frustration in z allows the system to form novel phases. 
Also frustration is expected to grow with increasing val- 
ues of the hopping amplitudes since the exchange cou- 
pling arises from the itinerancy of atoms. In the follow- 
ing, we shall address these two aspects by mapping out 
the ground state phase diagram for total filling one per 
site via BDMFT. 

Before going into a detailed discussion of the system's 
properties, let us at this point briefly introduce Bosonic 
Dynamical Mean Field Theory (BDMFT) 0-[|, which is 
an extension of Dynamical Mean Field Theory (DMFT) 
EH [3 ; originally developed to treat strongly correlated 
fermionic systems. It is non-perturbative, captures lo- 
cal quantum fluctuations exactly and becomes exact in 
the infinite-dimensional limit. Note that for antiferro- 
magnetic exchange J z the ground state of the system 
may break the translational symmetry of the lattice. 
We investigate this system using real-space BDMFT (R- 
BDMFT) p3, which is a generalization of BDMFT to 
a position-dependent self-energy and captures inhomo- 
geneous quantum phases. Within BDMFT/R-BDMFT, 
the physics on each lattice site is determined from a local 
effective action which is obtained by integrating out all 
the other degrees of freedom in the lattice model, exclud- 
ing the lattice site considered. The local effective action 
is then represented by an Anderson impurity model 0- 
.9]. We use Exact Diagonalization (ED) (TE E3 01 the 



effective Anderson Hamiltonian with a finite number of 
bath orbitals to solve the local action (nbath = 4 bath 
orbitals are chosen in the current work). Details of the 
R-BDMFT method have been published previously fl5j . 



III. RESULTS 

Since we are mainly interested in the effects of ge- 
ometric frustration, we focus on the parameter regime 
V a = Vb 3> U, where the leading exchange couplings be- 
come Ising antiferromagnetic for highly asymmetric hop- 
ping, resulting in Ising-type frustration. In our simula- 
tions, the chemical potentials fi a and \ib are tuned to 
equal particle number of both species (N a = X),:( n 'ia) = 
N b = J2i( n ib)) and a total filling factor p = Y^i( n ia + 
nib)/N\ at = 1 (iVi a t is the number of lattice sites.). The 
full range of interactions from strong coupling, deep in 
the Mott phase, all the way to the superfiuid at weak 
coupling is investigated within BDMFT. 

Our main results are summarized in Fig. [TJi, which 
shows the ground state phase diagram for large intra- 
species interaction strengths {V a ,b/U = 48). Two dif- 
ferent magnetic phases are found in the Mott insulator. 
When t a and tb are of comparable magnitude, the lead- 
ing exchange coupling is ferromagnetic in the icy-plane, 
and the system is in the xy-ferromagnetic phase, char- 
acterized by a uniform magnetization in the xy-plane. 
On the other hand, for sufficiently large asymmetry be- 
tween the two hopping amplitudes, we observe two types 
of spin-density wave (SDW) phases (see insets (a.l) and 
(a. 2) in Fig. [TJl) , which break the translation symmetry 
of the lattice. All of these spin-ordered phases are found 
to persist up to the superfiuid transition. For large hop- 
ping the ground state breaks U(l) symmetry and devel- 
ops superfiuid order (a),(b). Depending on the relative 
magnitude of t a and tb , the system can also exhibit addi- 
tional charge-density wave order in each species, leading 
to a supersolid. In the following subsections we discuss 
detailed properties of the different phases. 



A. Spin-density wave (SDW) in the Mott-insulator 
region 

Let us now discuss the origin of the two different SDW 
phases found in the asymmetric hopping regime. 



1. "3-sublattice" SDW 

In the less asymmetric hopping regime of the SDW 
phase region (yellow area in the phase diagram Fig. [T^) , 
the SDW features a 3-sublattice structure (see Fig. [5Jd) 
ordering at the wave vectors ±Q = (±47r/3,0). The 
emergence of this type of SDW can be understood from a 
spin- wave theory of the spin-I/2 XX Z model ©, which 
shows that the uniform a;y-ferromagnetic phase develops 
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FIG. 1: (Color online) (a) Ground state phase diagram at 
unit filling p = 1, with V a = Vb = 48(7, obtained from cal- 
culations on a 6 x 6 triangular lattice. We observe four ma- 
jor different phases, which are the homogeneous superfluid 
(HSF), ly-ferromagnet, spin-density wave (SDW), and super- 
solid (SS). In the SDW region, the area where the SDW has 
a 3-sublattice structure is marked in yellow, while the area 
featuring the stripe SDW is marked red. The insets (a.l) 
and (a.2) show the corresponding z-magnetization distribu- 
tion (Si). The small blue area indicates a coexistence region 
of HSF and stripe SDW phases, (b) The triangular lattice 
considered here with sublattice indices A,B,C marked near 
the corresponding lattice sites, ai and a2 are the two lat- 
tice vectors, (c) tj, dependence of the 3-sublattice SDW order 
parameter at fixed 9t a /U = 0.6. 

an instability at ±Q = (±47r/3, 0) towards 3-sublattice 
ordering with increasing J z /J± (la. fli^. Here, we choose 
the order parameter of the 3-sublattice SDW phase as 

4> = (S Z A ) + (S z B }e t2n/3 + {S z c )e l ^ /Z (6) 

where {S Z A B c ) is the z-magnetization on the sites of 
the A, B, and C sublattice respectively. 4> is the Fourier 
transform of the z-magnetization distribution (Sf) at the 
wave vector Q. In Fig.[Tp, at fixed t a , the tb dependence 
of \(f>\ is shown, indicating a second order phase transition 
from the 3-sublattice SDW to an xy-ferromagnet. 

Let us mention that previous investigations of the spin- 
1/2 XXZ model in the large J z /J± region have revealed 



that the z-magnetization favors a different 3-sublattice 
patten of ((S A ) = ±2m, (Sp ) = =F™, (Sq) = Tm) in 
the thermodynamic limit |2Ct l2l| , where m is a positive 
number characterizing the strength of the magnetization. 
However, although the effective spin-1/2 XXZ model is 
a reasonable description of the original two-component 
Bose-Hubbard model (fTJ) in the strong coupling limit, 
one can not exclude the influence of higher order terms 
neglected in it. As a matter of fact, even within the 
spin-1/2 XXZ model itself, although (±2m, =pm, =Fm) is 
favored in the thermodynamic limit, a metastable pat- 
tern (±m, Tfij 0) is a l s ° found in quantum Monte Carlo 
simulations on finite-sized lattices (20l . [2~j| . moreover a 
variational study shows that the energies corresponding 
to these two patterns are very close to each other (2^ . 
In our simulations, we find that the system is close to a 
(±m, =pm,0) configuration (see Fig. [2jj) , which could in- 
dicate that the effective spin-1/2 XXZ model plus higher- 
order corrections favors this pattern. 



2. Stripe SDW 

In the extremely asymmetric hopping region (area 
marked in red in the phase diagram Fig. [TJi), we observe 
that another type of SDW phase arises in the system. 
The z-magnetization distribution is then characterized 
by a stripe pattern, shown in Fig. [^F- The transition 
from the stripe SDW phase to the homogeneous super- 
fluid is of first order, which leads to small coexistence 
regions (blue areas in Fig.[T^). 

To understand the appearance of this stripe SDW, we 
notice that in this regime the hopping amplitude of one 
species dominates. Without loss of generality, in the fol- 
lowing discussion we assume t a ;g> tb, which indicates 
that the b species can be considered as almost immo- 
bile scattering centers. If we assume that the effective 
spin-1/2 XXZ model @ description still holds true in 
this case (in principle the neglected higher order terms 
may be relevant), the longitudinal exchange coupling J z 
is then much larger than the transverse one J±, hence 
the spin-1/2 XXZ model approximately reduces to an 
antiferromagnetic Ising model. From investigations of 
the antiferromagnetic Ising model on a triangular lattice 
[l[ , we know that the stripe pattern of the magnetization 
is one of an infinite number of degenerate ground state 
configurations. However, here we observe no other pat- 
terns in our simulation except the stripe configuration. 
This indicates that for quantitative insight one needs to 
go beyond the spin-1/2 XXZ model description, which 
we will do in the following. 

We notice that for extremely asymmetric hopping, to- 
gether with the condition that the intra-species interac- 
tions are much larger than the inter-species interactions, 
i.e. V a .b ^ U, the two-component Bose-Hubbard model 
([T|) can be simplified to a bosonic Falicov-Kimball model 
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with itinerant hard-core bosons 

H BFK = — 'y](t a alaj+h..c.)+Uy K /i riainbi— ^ p a n ai 

(i.j) i i;a—a,b 

(7) 

At zero temperature, integrating out the itinerant 
bosonic degree of freedom under the assumption t a <C U , 



we end up with an effective classical spin Hamiltonian 
representing the density-density interactions of the im- 
mobile b species, which is accurate to the order 0{t 5 a /U A ) 
(23| and reads 
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(8) 



where N\ at is the number of lattice sites, and 
denote next- nearest neighbor and next- next- 
nearest neighbor sites respectively (see Fig. (2ji). The 

Ising pseudo spin s, is defined as s, = (— 1)"> +1 . P de- 
notes the plaquettes made up of two triangles sharing an 
edge and sp = spiSp2Sp3Spi is the product of the spins 
on all sites in the plaquette. Detailed studies of the model 
© in Ref. [H] show that at the filling N a = N b = N/2, 
the ground state of the above classical spin Hamiltonian 
H®ff K (s; fi a ) is non-degenerate and has a stripe pattern 
configuration of the Ising pseudo spins similar to that 
observed in our numerical BDMFT simulations. From 
([5]) we observe that in the leading order 0(t^/U) the 
classical spin model shows an Ising-type frustration on 
the triangular lattice, however the frustration is removed 
by next-nearest neighbor and next-next-nearest neighbor 
effective spin interactions which originate from higher- 
order density fluctuations of the mobile particles (the a 
species atoms in this case). 



3. Weak charge- density wave (CDW) on top of SDW phase 

In the the SDW region we also observe a weak 
charge-density wave (CDW) of the total density pi 



,0.5 



<£o= 



i) forming on top of the SDW (see Fig. ^p). 



Specifically, we found that for t a > tb, those lattice sites 
with a negative ^-magnetization have a larger total den- 
sity pi, while for t a < tb, those sites with a positive z- 
magnetization have the larger density. 

The origin of this weak CDW can be easily understood 
by investigating a simplified two-site two-component 
Bose-Hubbard model 



H TT = -{t a a\a R + t b b [ I b R + h..c.) + U ^ 



-- v cxni a (n ia - 1) . 



i=L,R 



(9) 
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FIG. 2: (Color online) (a-c) z-magnetization distribution 
(Sf) for different Mott-insulating phases with quantum mag- 
netic order, (a) a;j/-ferromagnet for 9t a /U = 0.6 and Qtb/U = 
0.14. (b) 3-sublattice SDW for 9t a /U = 0.7 and 9t b /U = 0.1. 
(c) Stripe SDW for 9t a /U = 0.7 and 9t b /U = 0.05. In 
all the three cases, we set V a ,b/U = 48. (d) Schematic fig- 
ure where nearest neighbor, next-nearest neighbor, and next- 
next-nearest neighbor sites in a triangular lattice are indicated 
by the purple solid line, red dashed line, and blue dash-dotted 
line respectively. 



Since the onsite interaction is much larger than the ki- 
netic energy, we treat the hopping terms as perturba- 
tions. Since we are investigating the CDW on top of 
the SDW, we assume a simple symmetry-broken state 

^G^j = d L b R |0) as the unperturbed ground state, 

which implies that the L and R sites have equal total 
density (p^ = pp = 1) but opposite ^-magnetization 
((Sl) = ~(Sr) — 1/2)- If we take into account the hop- 
ping terms to first order, the ground state has the form 



i—L,R;a—a,b 



Ua{b R |0) + t b a\b\ |0) + t a a R b R |0) 

Vu 2 + tl + ti 



(10) 
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and the total density on the site L and R is 

PL « 1-6, (11) 
PR. ~ 1 + 6, (12) 

where 6 = (t 2 - tl)/(U 2 + t\ + t 2 .), indicating that the 
weak CDW order on top of the SDW originates from the 
asymmetry of the hopping amplitudes, i.e. the "lighter" 
species can more easily delocalize to neighboring sites. 
Fig. [3] shows the amplitudes of the stripe SDW and the 
CDW orders as a function of t a at fixed Those ampli- 
tudes are defined by dp = p+ — p- and S* tripc = — 
for CDW and SDW order respectively, where p_|_(p_) is 
the total particle density per site on the stripe with higher 
(lower) density and S+(St) is the z-magnetization per 
site on the stripe with positive (negative) value. 

More generally we can analyze the emergence of the 
CDW within a Ginzburg-Landau framework. We assume 
spin density wave order at wavevectors Q and — Q, and 
study its interaction with a putative CDW order param- 
eter. The lattice translation and Ising spin symmetries of 
the problem determine the allowed coupling terms that 
affect the CDW order, given in the following expansion 
of the free energy 

F = &s\Sq\ 2 + a pl \p Q \ 2 + a p2 \p2Ci\ 2 
+/3x(S^p_ 2Q +h.c.) 
+/? 2 |^| 2 | PQ | 2 + ft [(S1 Q ) 2 ( PQ ) 2 + h.c] 

+7s|Sq| 4 + 7piIpq| 4 +7 P 2|P2q| 4 , (13) 

where the a's, /3's, and j's are GL coefficients. The cru- 
cial point is that symmetry allows a linear coupling of a 
charge density wave at ±2Q to the square of the SDW 
order parameter (coupling (3\ above). Therefore a SDW 
order at ±Q will necessarily produce CDW at ^2Q- The 
opposite is not true, that is CDW order will not produce 
SDW in general. This is because a quadratic-linear cou- 
pling to the SDW is prohibited by the Ising spin symme- 
try. 

Let us first describe the implications of this fact in 
the case of the 3-sublattice SDW with wavevector Q = 
(±47r/3,0). In this case the CDW induced by the cubic 
coupling is of wavevector 2Q = — sgn(Q :r )(47r, 0) + Q, 
which in the lattice is equivalent to Q. This explains, 
from a very general argument, why a CDW at Q must 
be produced in this case. 

Let us move on to consider the case of the stripe SDW 
having a wavevector Q = ±K/2 with K being a recip- 
rocal lattice vector. The coupling to p 2 Q is then equiv- 
alent on the lattice to a coupling to the average density 
Pq=o that can be absorbed into the chemical potential 
term. As a consequence, the only non-trivial CDWs are 
the ones at ±Q and induced by the SDW order at ±Q 
through the next order terms quadratic in the CDW. 

To briefly summarize, both 3-sublattice SDW and 
stripe SDW induce CDW at the same wave vector, which 
is consistent with our simulations (see Fig. [3Jd, which 
shows a stripe SDW with a large SDW order induces a 




(a) (b) 

FIG. 3: (Color online) (a) In the stripe SDW phase region, we 
show the CDW order (blue solid line) and SDW order (black 
dashed line) as a function of t a when 9tt/U — 0.01 is kept 
fixed, (b) The total density (pi) (upper panel) and the z- 
magnetization ({Sf}) distribution (lower panel) for V a ,t/U = 
48, 9t a /U = 0.7, and 9t b /U = 0.01 . 



CDW at the same wave vector). However, in the case of 
a stripe SDW the CDW is induced by coupling to |pq| 2 
rather than by linear coupling. It therefore requires a 
critical strength of the stripe SDW to change the sign of 
the coefficient of \pq\ 2 and induce CDW order. 



B. Superfluid and supersolid in the large- hopping 
region 

At sufficiently large hopping amplitudes, the ground 
state of the system breaks U(l) symmetry and devel- 
ops superfluid long-range order characterized by non- 
zero values of (a) or (b). In the region where the hop- 
ping amplitudes t a and tb are comparable to each other, 
we observe a second order phase transition from xy- 
ferromagnetism to the homogeneous superfluid (HSF). 
On the other hand, for highly asymmetric t a and tb, in- 
creasing the hopping amplitudes can first drive the sys- 
tem from the SDW into a supersolid, which is charac- 
terized by coexisting superfluid and CDW order. In the 
supersolid phase, the lighter species develops superfluid 
order, i.e. (a) ^ 0, while the heavier species remains 
insulating, i.e. (a) = 0, where a(a) denote the annihi- 
lation operator of the lighter/heavier species and both 
species have CDW order in their density distribution re- 
spectively. When the hopping amplitudes are further in- 
creased, a transition from the supersolid to a homoge- 
neous superfluid is observed. 



IV. CONCLUSION 

We have investigated zero temperature quantum 
phases of Bose-Bose mixtures in a triangular lattice us- 
ing real-space Bosonic Dynamical Mean Field Theory. 
A rich phase diagram including xy-ferromagnet, spin- 
density wave, superfluid, and supersolid phases is found. 
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In the strong coupling regime, although an effective spin- 
1/2 XX Z model gives qualitative insight, interesting 
phases beyond this effective description are found: A 
stripe spin-density wave is identified for highly asymmet- 
ric hopping, which originates from the interplay between 
classical geometric frustration and higher order density 
fluctuations of the lighter species. Moreover, on top of 
the spin-density wave, due to asymmetric hopping am- 
plitudes, the system shows a weak charge-density wave 
in the total particle density distribution. 
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